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00 . We consider the dynamics of a system coupled to a thermal bath, going 

o, 

beyond the standard two-level system through the addition of an energy exci- 
tation degree of freedom. Further extensions are to systems containing many 
fermions, with the master equations modified to take Fermi-Dirac statistics 

into account, and to potentials with a time-dependent bias that induce reso- 
rt , 

nant avoided crossing transitions. The limit Q — > oo, where the interaction 
rate with the bath is much greater than all free oscillation rates for the sys- 
tem, is interrogated. Two behaviors are possible: freezing (quantum Zeno 
effect) or synchronization (motional narrowing). We clarify the conditions 
that give rise to each possibility, making an explicit connection with quantum 
measurement theory. We compare the evolution of quantal coherence for the 
two cases as a function of Q, noting that full coherence is restored as Q — > oo. 
Using an extended master equation, the effect of system-bath interactions 
on entanglement in bipartite system states is computed. In particular, we 
show that the sychronization case sees bipartite system entanglement fully 
preserved in the large Q limit. 
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I. INTRODUCTION 



The dynamics of a quantum system coupled to a thermal bath is of great interest in 
many branches of physics [|I| . This topic is of practical importance because no system apart 
from the universe as a whole is actually isolated. The contemporary focus on quantum 
computation |2| has been providing added impetus to these studies, as has the need to 
understand collision-affected neutrino oscillation dynamics in astrophysical and cosmological 
settings ||. From a fundamental perspective, the effect of system-bath interactions on 
quantal coherence in the system bears on the measurement problem and on the emergence 
of classical behavior out of an underlying quantum dynamics @]. 

In this paper we will revisit two closely related phenomena that arise in the limit of very 
rapid system-bath interactions: freezing (or quantum Zeno effect) and motional narrowing 
(or synchronization). Our intentions are twofold: (i) To present a unified perspective on these 
effects that emphasizes their connection with each other, and their connections to general 
quantal coherence and measurement questions, (ii) To extend the analysis to incorporate 
Fermi-Dirac statistics, a potentials with time dependent bias, and 2-particle entangles states. 
For definiteness we will focus on a double-well system, which turns out to be a versatile 
container for exemplifying and coordinating a number of phenomena. However almost all of 
our considerations can be carried over to other arenas, e. g. spin systems. 

A. Freezing or Zeno 

What we call freezing in this paper, the drastic slowing of a process through environ- 
mental interactions, can be discussed in a number of ways. If we start with an oscillator of 
some kind, and crank up a damping parameter to the point that the oscillator hardly budges 
in any reasonable amount of time, we might call the system overdamped. If we start with 
a decaying system on which we make very rapid measurements to see whether the original 
state is still there and find that the system has hardly evolved, we call it the quantum Zeno 
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effect (QZE) [||]. Except for the fact that something has been slowed, two behaviors that 
are described so differently appear to have little in common. But when the "over-damped" 
behavior in the oscillator comes from interactions with a quantum thermal bath, and the 
"observation" in the QZE is replaced by interactions of the system with such a bath as 
well, then we find that the two terminologies can describe the same essential physics. If we 
particularize the "oscillator" in the above to the case of a particle in a symmetrical double- 
well, then it turns out that we must require the bath interaction to distinguish left from 
right in order for freezing to occur; hereafter we shall refer to such a bath as "side-sensing" . 
Correspondingly, for the particle-decay prototype of the QZE, in order to simulate rapid 
measurements the bath must sense the presence or absence of the initial particle. For the 
widely discussed two-state systems, the freezing phenomenon is essentially the same whether 
they are double- well systems originating in a condensed matter context |j] , molecules || [|7| , 
nuclear spins, squids ||, or neutrinos in the early universe 0. 

In practical circumstances we shall want to know the answer to quantitative questions 
such as how frozen a system becomes as a result of bath interactions, the answers to which 
can be roughly expressed in terms of ratios of the different time scales in the problem. We 
also want to know how the slowing or freezing behavior is correlated to other calculable 
phenomena in the system. 

B. "Synchronization in a linear system" or "motional narrowing" 

These two effects seem to be identical at the formal level. They carry different names as 
a result of arising in very different physical contexts. Synchronization is a natural partner 
to the slowing or freezing phenomena for the case in which the bath does not distinguish 
which of the two states the two-state system is in. We shall describe such a coupling as 
"side-blind". For a pure two state system, this describes an interaction that does nothing 
at all of note. But in our example of an electron in a symmetrical double well we may 
wish to place the system in a bath that can excite or deexcite the states that are higher 
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lying than the first two (i.e. higher lying than the nearly degenerate symmetrical and anti- 
symmetrical combinations of the single- well ground states). We shall refer to the inclusion 
of these states as the introduction of a "vertical" coordinate, using "horizontal" to denote 
the left and right coordinate. The oscillation rates are different for the higher-lying states, 
so that the density matrix for a particle that is, say, initially on one side of the well no 
longer oscillates sinusoidally between the two states. But in the limit of large "side-blind" 
coupling to the bath the simple sinusoidal behavior is restored, as though the different 
frequencies had been made the same, and the oscillations of the different vertical components 



synchronized. Similar behavior in a neutrino-based example has been reported in ref. [10 



"Motional narrowing" has been used to describe analogous behavior in an NMR system 
|TT| . Qualitatively, motion narrowing occurs when the system is made to change state so 
rapidly that it is only able to respond to an average Hamiltonian. That is, the "horizontal" 
evolution is governed by parameters which are averaged over the "vertical" states. 

In this paper we review the basic analysis required to establish the above assertions, 
including the derivation of a master equation, all in the double- well context, and we provide 
some plots of typical behavior of the solutions. Then we investigate a number of new topics 
related to the limits mentioned above: 



1. Several fermions in the well 

This entails writing a modified master equation to take into account Fermi statistics. 
The solutions, in the limit of large "side-blind" noise, have a counterintuitive property: an 
assemblage of several particles can move sinusoidally from one side to the other, giving the 
appearance of being in the same state. 



2. Systems with four horizontal states 



This is like the familiar case of two spin degrees of freedom; in our well realizations we 
utilize two independent particles in the well, each particle with the same structure of vertical 
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levels. We can now study entanglement of the horizontal states. In the presence of moderate 
"side-blind" noise we find rapid destruction of entanglement, while in the "synchronizing" 
limit of large "side-blind" noise we have, in contrast, perfect preservation of entanglement. 
The latter may find application in quantum information studies, just as its freezing cousin 
has already done |T2| ]. 

3. Avoided crossings 

If we add a variable bias to the symmetrical well, so that one side is deeper than the 
other, we can perform interesting numerical experiments in which the bias is slowly changed 
from one sign to the other. In the absence of the coupling to the thermal bath, a particle 
localized on one side can be permanently transported to the other side through an adiabatic 
reversal of sign of the bias. When moderate "side-blind" noise is added, we find that instead 
of the near 100% transportation of the uncoupled case we obtain a final distribution with 
approximately 50% of the probability on either side. However, in the presence of strong 
"side-blind" noise we find the efficiency of the transportation to be restored. 

4- Quantum measurement 

We discuss the connection between system-bath interactions and quantum measurement. 
We show that the system-apparatus-environment entanglement needed within the emerging 
paradigm for the description of quantum measurements can arise naturally in our context 
for the side-sensing case. 

C. Structure of the paper 

In section |T| we specify our system-bath Hamiltonian and the corresponding master 
equation for the system. We then discuss both the synchronized and frozen behavior and 
include an extension to the case of multi-particle systems with Fermi statistics. At the end 
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of this section we analyse the case of the double well with a time dependent bias, focussing 
on avoided level crossing transitions. In section |T| we determine the evolution of entangled 
well states under the influence of side-blind noise. We discuss the side-sensing and side- 
blind interactions from a measurement theoretic perspective in section [TV] and in section [V] 
we conclude. 

II. DYNAMICS OF FREEZING VERSUS SYNCHRONIZATION 

A. The master equation: single particle case 

We focus on those situations where the effect of system-bath interactions on the dynamics 
of the system can be derived from a master equation for the reduced density matrix of the 
system. 

The horizontal "spin" coordinate of the n-state system (where n = 2 for the double well) 
is indexed with Greek letters, while the vertical energy variables will be denoted E^. The 
system basis states are denoted \E i: a). The dynamics of the system uncoupled to the bath 
is governed by the Hamiltonian, 

tf sys = W«p + A a/3 (^)]|^,a)(^,/3|. (1) 

i,a,/3 

We will study the case \ a p{Ei) -C Ek for all values i, k. This means that the energy 
eigenstates fall into nearly degenerate pairs near each Ei, and the energy splittings within 
these pairs are always much smaller than the spacings AE between the pairs (single-well 
spacings) . 

Consider a reduced density matrix for the system p sys (t) = Trb a th[p s+b (^)], where p s+h (t) 
is the complete density matrix of the model. At an initial time t = we choose a form 
diagonal in the vertical indices, 

p ^(t = 0)=J2 Pa f s(E l ,t = ti)\E h a)(E h {3\. (2) 



6 



For these initial conditions, one of the results of a master equation derivation is that for times 
t ^> (AE)^ 1 , the operator p sys remains so nearly diagonal in the vertical (energy) indices that 
we can continue to describe the system by a vector in the energy space, p(Ei, t) a /3, replacing 
(t = 0) by (t) in Eq.@. Henceforth we suppress the indices (a,/?) in the horizontal space, 
in which p{E il t) remains an operator. 

We now have to supply a model for the system-bath interaction. A natural and simple 
case to consider is a factorised form, 

H s - b = (V, (3) 

where ( operates on the horizontal indices and V depends only on the vertical coordinates 
and the coordinates of the bath. This coupling is also taken to be weak, in a way that allows 
the problem to be solved, although it can at the same time be strong in another sense: we 
will eventually study the limiting case where the system-bath interaction rate is much larger 
than the free oscillation rates of the system. 

In the Appendix, we show that for this form of coupling the master equations, generalized 
to include the vertical structure, are of the form, 

^p(E h t) = -i[\(Ei), p(Ei, t)]+J2( P(Ej, t) Cr(%, Ei) 

- \ (C 2 p(E h t) + p(E h t) C 2 ) £ T(E h Ei). (4) 

Z 3 

These equations are of Bloch form. All of the elements in Eq.(^) are matrices in the horizontal 
space except for the rate functions V . They obey the relation 

Y{E h Ei) = exp[(^ - ^)/T]r(^, E d ) (5) 

due to the thermal equilibrium of the bath. For the neutrino application, and with the 
right identifications, Eq.(f|) gives the "quantum kinetic equations" derived by McKellar and 
Thomson |Tj|. They can also be found in the NMR literature, for example in Ref. P| . The 
derivation in the Appendix produces the equations in the context needed for our present 
results. 
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We now specialize to the symmetrical double well. The two-dimensional horizontal space 
consists of the left and right sides of the well, conveniently labelled by the eigenvalues of cr 3 , 
namely +1 (—1) for the particle to be on the left (right) side. The barrier height is chosen 
to be large compared to AE. The vertical states are the single well energy excitations. We 
take 



in order to model the free oscillations between left and right. The energy splittings g(Ei), 
which are essentially tunneling rates, have strong Ei dependence. For the "spin" dependence 
of the coupling to the bath we take, 



where the parameter 6 plays the important role of modelling the horizontal structure of 
the system-bath interaction. In the case 6 = 0, the transition matrix elements of the 
operator H s _ h for Ei — > Ej in the left-hand well, considered by themselves, are the same 
as those for the right-hand well, as would be the case in the dipole approximation of the 
interaction with a radiation field. For the case 6=1 (—1), only the amplitude on the 
left-hand (right-hand) side of the well interacts with the bath. In other words, 6 = (±1) 
describes left-right symmetric (maximally asymmetric) system-bath couplings: side-blind 
(maximally side-sensing) in other words. As we review below, frozen behavior arises as a 
limiting case for 6^0 couplings, while synchronization or motional narrowing occurs as a 
limiting case when 6 = 0. 

The function T(E i) Ej) is the (single well) rate for a state Ei to make a transition to state 
Ej. When T ps AE, the bath induced transitions among the vertical states are dynamically 
important. In this case our models differ substantively from models that do not have the 
vertical structure. 

As an example we use a double square well with infinite barriers to the left and right 
of x = ±8a, and with a central barrier of width 2a and height Uq. We choose Uq and the 



Kp(Ei) = [a^\ a pg{Ei) 



(6) 



C = 1 + ba 3 , 



(7) 
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particle mass such that twenty states are bound with energies less than Uq. The low-lying 
states are thus very nearly the symmetrical and anti-symmetrical combinations of separate 
states in an infinite well of width 7a. The energy splittings, g(Ei) are easy to calculate in 
the limit of small tunneling. 

For the bath we take bosons with uo = vk, equal energy spacing 2nv/L where L — > oo, 
and with creation and annihilation operators, a/ n , a n . We take 

H = J2^ n aia n + H^ s \ (8) 

n 

and in the system-bath interaction, H s ^ = V(, we take the vertical operator V to have 
off-diagonal matrix elements given in dipole expansion as 

V ik = 7 HK a n + a n)5 jk + ik n x jk (a ] n - a n )\. (9) 

where q is a coupling strength. We calculate the off-diagonal dipole matrix elements, Xjk, 
using the infinite well wave- functions and define Cjk = qk n Xjk, where the later formalism 
will justify using (near) energy conservation to express k n as ±(-Ej — E^)/v. The transition 
rates in the bath interaction part of Eq.@ are then, 

, , r d(^Ej — Ek) 0(Ek — Ej) "1 . . 

T{E k , Ej) = 2nc jk c kj [ e(E ._ Ek)/T _ - + - _ e{Ej _ Ek)/T \ ■ (10) 

We will now display solutions to the master equations for the cases 6 = and b ^ 0. 

B. Results for 6 = 0: synchronization or motional narrowing 

We first consider side-blind system-bath couplings, so that 6 = 0, and we solve Eq.(f|) 
for several values of q 2 . The resulting curves are labeled with a ratio Q that is proportional 

to q 2 , 

Q = (r)/(g), (ii) 

where (g) is the thermal average of the oscillation rates, g(E), and (r) that of the bath 
induced transition rates, T(Ej, E k ). Figure 1 shows the probability that a particle beginning 



on the left side in a thermal distribution will be on the left side at time t. We choose 
temperature T = 5 in units of the single-well ground state energy. 




500 1000 1500 2000 2500 3000 

FIG. 1. The probability P for a particle to be found in the left hand well, for a bath temperature 
T = 5, in units of the ground state energy of the single infinite well. The initial condition is that 
the particle is on the left with a thermal distribution of energies. The curves are labelled by values 
of the rate-ratio parameter Q that range from 0.001 to 1000. The time is in an arbitrary unit, the 
scale of which is of the order of fi times the inverse of the single well energy spacings. 

For the largest value of Q the probability shows undamped sinusoidal behavior in time. 
(This can be demonstrated in the Q — > oo limit analytically, with the result that the oscil- 
lation frequency is the thermal average of the oscillation frequencies for the separate modes 
||10||). For somewhat smaller values the probability settles rather quickly to 50%. For very 
small Q the probability shows the irregular behavior characteristic of the incommensurate 
frequencies g(Ej). Note that arbitrarily large values of Q can be achieved for any value 
of system-bath coupling, q (whether the mechanism be photons, phonons or collisions with 
surrounding particles), by widening sufficiently the barrier between the wells. 

In Fig. 2 we show the von Neumann entropy, as defined by 

S = -^{p(E s )]og\p(E j )]} 1 (12) 

3 

where the trace is over the two-dimensional horizontal space. We note that as Q is increased 
from small to moderate values the rate of entropy increase rises, but that as Q becomes very 
large the rate of entropy increase goes to zero. 
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FIG. 2. The entropy S as a function of time for the same values of Q as used in fig.l. Time in 
arbitrary units. 

Notice that the entropy begins at a nonzero value. This is simply due to our initial 
state being thermal with respect to the vertical variable. (If we had chosen a non-thermal 
state it would have been rapidly thermalized anyway in the large Q limit.) Any increase in 
S beyond this base level then reflects loss of order with respect to the horizontal variable, 
which is really what we are interested in. 



We briefly summarize some examples. In the case plotted in Fig. 1, for the strongest 
coupling, we find little change for 6 = 0.001. But 6 = 0.005 starts to poison the oscillation, 
damping it by a factor of 20% in five periods. As we move to much larger b, the motion 
becomes more frozen, so that when, for example, b = 0.5, 98% of the probability has 
remained on the left hand side over a time span of five average periods of the 6 = system. 
Thus side-sensing bath interactions are very efficient in freezing the time development. 

The von Neumann entropy S shows qualitatively similar behavior to the 6 = case. As 
we increase R, the entropy at first increases, but eventually it becomes frozen at its initial 
value when the QZE limit is reached. An important difference between the two cases is that 
entropy increase abates for 6 = with the system evolving non-trivially with time, whereas 



C. Results for 6/0: freezing or quantum Zeno effect 
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for b 7^ it freezes, simply because all evolution freezes. 



D. Fermi statistics case 

We presented the results for the previous sections in terms of the density matrix for a 
single particle in the well. If we choose to put several distinguisable particles in the same well 
and there are no interactions among these particles, then the results for each particle are the 
same as the previous results, and if we wished we could write a single equation for the density 
matrix for the assemblage, normalized to the number of particles. The same is true for the 
case of an assemblage of indistinguishable particles in the case that the temperature is high 
enough for Boltzmann statistics to obtain. For the iV particle case with Fermi statistics, 
in the regime where the statistics matter, we have not succeeded in deriving an acceptable 
modification of Eq.(^). However, if we choose to use a formulation in which the diagonal 
matrix element of a density operator gives the probability that a state is occupied, and these 
occupancies are taken as independent [], there is a simple extension of Eq. that describes 



the time evolution of the ensemble of fermions in the well. Following the argument of |n 
we replace Eq.(H) by, 

&p(E i ,t) = -i[\(E i ),p(E i ,t)] + 
\ Ej([l - p(Ei, *)]Cp(#i, *) C + C p(Ej,t) C [1 - p(E, t)])T{Ej, Ei) 
- \ Ej (C[l - p(Ej, t)}( p{E h t) + p(E h t) C[l - p(E d , t)]C)T(Ei, E 3 ). (13) 

The density matrix is now to be normalized to the average number of particles, 
E, T^lftiEj, t)] = (N), which is time independent. As an example we take a system which 
begins at t — with no particles on the RHS, and an average number of particles (N) = 8 
on the LHS, for a well system with the same parameters used to generate data for the single 



1 When the occupancies are taken as independent, the particle number is necessarily indeterminate. 
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particle calculation. For the initial energy distribution we take 50% occupancy for each of 
the lowest sixteen states. 

In the case of zero coupling to the thermal bath, the various particles rattle back and forth 
independently with their incommensurate frequencies. Our terminology has been changed 
here, since previously we pictured the density matrix for a single particle performing exactly 
the same gyrations. We note that in the case without vertical transitions the Fermi statistics 
enter only in the restrictions put on the initial state. 

In the presence of a moderate coupling to the bath, the system tends rather quickly to 
an equilibrium state with 50%-50% L-R occupancy, as shown in the heavy dashed curve in 
Fig. [|. In this state the final energy distribution is a Fermi function, 

p(E j ,t) = [l + e&-M T ]- 1 (14) 

where the parameter \x is set by 

2X;[l + e (Bi "" )/T ]" 1 = (N) (15) 

3 

But, as before, if we increase the bath coupling sufficiently, the particle density oscillates 
back and forth sinusoidally as shown in the heavy solid curve of Fig. [| The average fre- 
quency is now given by the average of the individual-state rates with respect to the Fermi 
distrubution function. We have not found in the literature an analogue of this behavior, 
in which a group of fermions undergoes collective synchronized oscillations, as though they 
were all in the same state. 
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FIG. 3. The time development of the average number of particles in the left hand well. The 
initial conditions are: the first sixteen levels of the left well each occupied with 50% probability. 
The heavy dashed curve gives the behavior with extremely small coupling to the bath, the lighter 
dashed curve the behavior with an intermediate value of the coupling, and the heavy solid curve 
the behavior with the largest coupling. The respective couplings are in the ratios of 1CT 3 , 1, and 
10 3 . The time is in arbitrary units. 



A generalization of the above considerations is the extension to a double well with a time 
dependent bias, 



Consider a case in which we begin with e > and e ^> \g(Ei)\ for some set of low-lying 
states (but e <C AE). A particle not interacting with the medium, initially localized on 
the left, with its density matrix distributed among the low-lying states, will stay largely on 
the left as long as the bias is maintained at the initial value. If the bias is slowly lowered, 
going through zero and ending at the negative of the original value, then the particle will 
efficiently be carried to the right hand well. The mechanism can be described as "avoided 



E. Adiabatic inversion 




e{t)a 3 . 



(16) 
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crossing" . [] Now we consider the same process in the presence of a side-blind thermal bath. 
With a moderate coupling to the bath, the system will still remain localized on the left if 
we keep the bias constant at its initial large value. But if we proceed as before through 
an adiabatic reversal of the bias, then the jostling caused by the vertical transitions will 
reduce the efficiency of the left to right transportation, and will tend to leave one-half of the 
probability behind in the left hand well. 

However, if we increase the strength of the bath coupling sufficiently, then we regain the 
efficient shift from left to right Figure |4] shows the results of the calculation, for the cases 
of no bath coupling, moderate bath coupling and strong bath coupling. Note that in this 
example the transport of the particle from the left to right is even more efficient in the case 
of large bath coupling than it was in the case of no bath coupling, the small inefficiency 
caused by some non-adiabaticity in the bias-changing process being reduced in the presence 
of the bath coupling. 

This rather surprising behavior can perhaps be understood qualitatively in the following 
way: First consider the previously discussed synchronized oscillation limit for the symmetri- 
cal well. There the synchronized frequency is somewhat higher than we would have guessed 
from the oscillation frequency for a typical particle (with E rs T ) in the absence of of the 
bath, because in our model the tunnelling rates increase rapidly as we move up in energy. 
It is also true that the non-adiabaticity, as defined by the fraction of particles ending up on 
the left at the end of our bias-reversal, declines rapidly with increasing energy, again in the 
absence of bath interactions. Thus the bath's effects, still a bit mysterious, are to use the 
synchronized frequency for the purpose of the estimation of adiabaticity. We emphasize that 



2 This is the analogue of the Mikheyev-Smirnov-Wolfenstein (MSW) |15| transition in neutrino 
physics. This behavior is addressed in || for the case of a thermal bath that does distinguish the 
right and left sides of the well. As expected, their results are quite different for the case of large Q. 

3 Detailed analysis of similar behavior for a neutrino application is given in [10]. 



15 



the real energy distribution of the particles remains very nearly thermal; but synchronization 
lets us do the tunnelling, in effect, at a higher energy, with consequent smaller inefficiency 
This phenomenon offers some promise for obtaining greater L-R transition efficiency in 
an application. Another example of parameters, beyond the ones presented in Fig. |] is the 
following: take the temperature so low that the occupancy of the ground state is 97%, the 
first excited state, 3% and the higher states negligible. Then we find, for a particular rate 
of bias inversion and a particular level of coupling to the bath, that there is 6% probability 
to end in the left hand well. This is to be compared to 16% in the absence of the bath 
coupling. 

p 



l 




5000 10000 15000 20000 25000 30000 35000 

FIG. 4. The probability that the particle is on the left, where the well bias is reversed. The 
heavy dashed curve is for zero system-bath coupling, the light dashed curve for moderate coupling 
and the solid curve for large coupling. The time is in arbitrary units. (The dependence of e on t 
has been taken to be proportional to (t — to) 3 , where to is the zero-bias time.) 

III. EVOLUTION OF ENTANGLEMENT 

Irrespective of the value of b, we have concluded that the quantum coherence of a state 
that is initially a left or right eigenstate (with a thermal distribution of energies) is preserved 
under time evolution in the Q — > oo limit. For b 7^ 0, this is because of freezing, while for 
b = coherence manifests as motional narrowing or synchronization. 
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It is natural to then ask whether entanglement is preserved in the same limit. For the 
QZE case, the answer is no, simply because any L-R entangled state has off-diagonal entries 
that are rapidly damped when & ^ 0. In other words, entanglement will be destroyed, 
not frozen, in the QZE limit. We now want to see what happens in the synchronization 
case. Intuitively, one expects preservation to be maintained because all components of the 
entangled multi-particle system will oscillate between the wells at the same rate.0 This 
intuition is borne out by direct computation. 

We will focus on bipartite or two-particle states for simplicity. The maximally entangled 
or Bell states are 

KT/± P P N = \Ej, L) ® \Ej, R) ± \E U R) ® \Ej, L) 
I * j &j) — ' 

|(f)± jp \ \Eu L) ® \Ej, L) ± \E t , R) g \Ej, R) 

|Q ,Ei,Ej) = -j= , (17) 

where one well particle has energy Ei and the other energy Ej. 

A straightforward generalisation of the derivation sketched in the Appendix yields the 
master equation for a general two-particle system. We denote the basis states of the sys- 
tem as \Ei, E'j,a, (3) and consider a density matrix (diagonal in the two energy spaces) 

Suppressing the flavor (R-L) matrix indices, the Bloch equations, generalized to include 
the vertical structure, are of the form, 

j t p(E h E 3 ,t) = —i[X(Ei),p(Ei, Ej, t)) — i[X'(Ej),p(Ei, Ej, t)) 

+ p(E k , Ej, t) (T(E k , E,^ + £ C' p(Ei, E m , t) ('T(E m , Ej) 

Ek E m 

- \ (C 2 p{E h Ej, t) + p{E h Ej, t) C 2 ) v ^ E k) 

1 ''^ 2 p(E i ,Ej,t)+p(E i ,Ej,t)((') 2 )Y,r(E j ,E m ). (18) 

E m 



4 Note that synchronization preserves coherence given any initial state, as opposed to the QZE 
which preserves coherence only if the initial state is an eigenstate of 03. 

17 



Here \(E), X'(E) , (, and (' are matrices in the (R-L, R-L) space. For the case of 
symmetric double wells, the only horizontal matrices in the free oscillation part of the 
Hamiltonian are constructed from the o\ matrices of the two particles and we have 

\' aa , ppl {E')=g{E')5 aal [a 1 }w, (19) 

It is useful to parametrise the the density matrix in terms of the coefficients x^ in an 
expansion in terms of sigma matrices 

p = x ij (E,E')a i ®a j . (20) 



If we choose as our initial condition one of the four Bell states in eq . fll7|) , we need consider 
only the six components Xoo, Xu, £22, £23, £32 and 233. If, additionally, the initial state is one 
of thermal equilibrium with respect to the vertical variables then xqq(E, E') and xn(E, E') 
remain constant and we are left with the four equations 

j.Xij(E, E') = g(E)Y,&kX kj (E, E') + g(E')J2tjkX ik (E, E') 

k k 

+ £ [x hj (E k , E' p t)T(E k , Ei) - Xij {E u E' p t)T(Ei, E k )_ 

+ [x ij (E i ,E' m ,t)r{E' n ,E' j ) -x^E'j^TiE^E'j] , (21) 

where the nonvanishing elements of the matrices £ are £32 = —£23 = 2. In eq.(pl]) we have 
specialised the collision terms to the case £ = £' = 1, in which the system-bath couplings do 
not distingiush the value of the (L,R) index. 

We quantify the degree of entanglement using the entanglement of formation ||16|| . For a 
bipartite system, the entanglement of formation is given by 

E f (C) = h(l+ 1 -VT^), (22) 

where 

h(x) = — x\og 2 x — (1 — x) log 2 (l — x). (23) 
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The "concurrence" C is 

C = max{0, Ai - A 2 - A 3 - A 4 }, (24) 

where the A, are the square roots of the eigenvalues of the matrix pp, and the spin-flipped 
state p is defined as 

p = (cr y ® a y )p*(o y ® a y ). (25) 

We illustrate the evolution of the entanglement in the figures below, for the side-blind 
case with the initial conditions |$ ± (£'i, E\)). In the absence of the heat bath, a system 
initially in the state \<&~(Ei, Ei)) will remain in that state. An initial state of |$ + (i?i, Ei)), 
however, evolves as \ij>(t)) = cos[2#(£i)t]|$+(£i, £1)) - i sin[2#(£i)£]|^+(£i, Ex)) in the 
absence of a heat bath, while, of course, maintaining maximal entanglement. 

In Figs. [| and[7] we plot the evolution of X 33 , 

X 33 = 2 ^0.25 + ^ x 33 (E, E 5 ) \ , (26) 

which measures the probability that if one of the particles is on the left (right), the other 
particle will also be found on the left (right). In Figs. |6] and [8] we plot the entanglement of 
formation, calculated after first summing over the vertical states 

P = J2p(E 1 ,E j ). (27) 

X33 




<g>t 



FIG. 5. The evolution of X33 as a function of time, for the initial state \& + (Ei, E±)). The 
dashed, dotted, solid and dot-dashed curves are for Q = 0.0005, 0.05, 5 and 500 respectively. 
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FIG. 6. The entanglement of formation, Ef, for the same parameters as FigJE]. 



X33 




<g>t 



FIG. 7. The evolution of X33 as a function of time, for the initial state |<J> (Ei,E\)). The 
values of the parameter Q are as in Fig.||. 



E f 
1 tr- 




<g>t 
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FIG. 8. The entanglement of formation, Ef, for the same parameters as in Fig.^. 

For small to moderate Q values, we observe the gradual relaxation and decay of entan- 
glement arising from the population of energy levels with different frequencies. For large 
values of Q, we obtain synchronised behaviour, where the evolution of the system resembles 
that of two isolated qubits, with X(Ei) replaced by an averaged frequency A. That is, despite 
the strong coupling to the bath, the system behaves, effectively, as though it were isolated. 
In the limit Q — > oo, we see the persistance of maximal entanglement. 

Finally, note that entanglement may be created or destroyed if the Hamiltonian were 
to have a term describing an interaction between the two qubits. In the large Q limit this 
process will also proceed as if the system were isolated, with the Hamiltonian parameters 
replaced by their averaged values. 



IV. DISCUSSION 



The freezing phenomenon or QZE is well known. It should be noted, though, that 
much of the recent non-neutrino literature does not study systems with a vertical structure. 
Indeed, when the vertical structure is removed Eq.(£|) reduces to the familiar form, 

dp . r , , r, 



Of 



-i[\,p}--[(, [C,p]}. 



(28) 



The decoherence term is then simply 



2 [6 



-2b 2 T 



( p 12 ^ 



(29) 



\P*12 J 

from which freezing can be deduced in the large Q limit, provided that b ^ 0. The presence 
of vertical states evidently does not alter these qualitative conclusions: freezing occurs at 
large Q for b ^ whether vertical structure exists or not. 

The above derivation also shows that 6 = system-bath interactions have no affect on 
the system in the absence of vertical structure. But, when it is not absent, we have seen that 
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there is an interesting conceptual cousin to the b ^ quantum Zeno effect: synchronization 
(or motional narrowing) when b = 0.0 Both occur in the Q — > oo limit, when the dynamics 
of the system is dominated by its interactions with the bath. While synchronization from 
motional narrowing is a known phenomenon, we have not found a clear and comprehensive 
discussion in the literature about its relationship with the QZE, quantal coherence and 
measurement, so we now provide one. 

The strong connection between freezing/synchronization and coherence is revealed by 
the evolution of entropy. Pure states always have vanishing entropy. Decoherence in the 
system induced by system-bath interactions can cause initially pure states to evolve into 
mixed states. In this case the entropy increases towards a maximum which is attained when 
all states in the statistical mixture are equally populated, with complete decoherence. 

We have to be careful to distinguish between what one may call vertical and horizontal 
contributions to entropy. As explained when we discussed Fig. 2, the entropy of our system 
starts out being nonzero because of the thermalization in the vertical variable. We are 
really interested in whether or not S increases from this initial offset, because that change 
would be caused by growing horizontal disorder. But we have seen that a system initially 
in a horizontal eigenstate experiences no entropy increase for any value of b in the large Q 
limit. This means that the system evolves in a pure-state manner (subtracting off the initial 
vertically induced offset), despite having its dynamics dominated by interactions with the 
bath: environmental decoherence hits a maximum at an intermediate value of Q, not in the 
limit Q — > oo. For b ^ 0, this is in a sense Pyrrhic victory over decoherence since it comes 
at the cost of freezing. For b = on the other hand, the ensemble described by the reduced 
density matrix evolves nontrivially as if it were a single isolated particle undergoing unitary 
evolution at an averaged frequency. 



Of course this manifests itself only when the vertical coordinate exists because one cannot talk 
about synchronization for a one-oscillator system! 
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Whether freezing or synchronization results in the large Q limit depends on how the 
bath couples to the system. If the bath is side-sensing, (b 7^ 0) then freezing happens; 
if it is side-blind (b = 0), then synchronization occurs. It has often been remarked that 
the QZE case is connected with the measurement process. Since the environment has a 
b 7^ coupling, it can "tell the difference" between left and right sides of the well, and thus 
"measure" this observable (the Hermitian operator 03). When Q — > 00, the system-bath 
"measurement" interactions happen repeatedly and quickly (constant monitoring), forcing 
the system to remain in the same 03 eigenstate (or more generally to remain in the same 
statistical mixture of eigenstates). On the other hand, if the coupling is of b = form, 
then there is no sense in which the environment probes whether the system is in the left or 
right state, and therefore it cannot "measure" and thus freeze 03. The oscillations between 
left and right are undamped, but the system states transform between vertical Ei levels so 
rapidly that they effectively all oscillate at the same averaged frequency. 

These are heuristic observations, as indicated by the scare quotes around key phrases and 
words. It would be useful to make the connection between b and measurement more rigorous. 
This should be possible in light of the recent clarification of how exactly measurement can 
be modeled in quantum theory. 

The emerging paradigm connects the measurement process to entanglement between sys- 
tem, apparatus and environment states [Q. We will see that microscopic collisions between 
a system particle and a spin-0 boson from the bath can be interpreted as inducing this sort 
of entanglement in the b 7^ case.Q So, the microscopic perspective will provide us with 
a simple explanation for why the b 7^ and 6 = cases are different, and indeed why the 
former is connected with measurement. These considerations will not capture the full com- 
plexity of the phenomena generated by the master equations (synchronization in particular), 



3 System-environment entanglement should not be confused with the intra-system entanglement 



studied in Sec. ID. 
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but they will furnish some physical insight. 

We first briefly summarize the emerging paradigm. The world is divided into three 
parts: system, apparatus and environment. Suppose that the system is initially in state 
\<S) = J2 n Cn(0)|n), where \n) is an eigenstate of O, the system observable that the apparatus 
can measure. Initially, we take the apparatus and environment to be in states |^4) and \£), 
respectievly. We also suppose that initially there are no correlations between the three 
subsystems, so the state of the world begins as simply |<S) <g> \A) <S> \S). By time t±, the 
system and apparatus have unitarily evolved into the entangled state J2 n c n(ti)\n) <E> \a n ), 
so that a correlation now exists between eigenstates of the system \n) and some basis states 
\a n ) for the apparatus. In the typical picture, the environment at this "pre-measurement" 
stage remains uncorrected with the system-plus-apparatus hybrid. The measurement is not 
complete at this stage because quantal correlations exist between the "outcomes" . By time 
t 2 , however, we suppose that unitary evolution has also entangled the environment with the 
apparatus, in a way such that each apparatus state \n) is linked to a distinct environmental 
state, \e n ). Now the world is in the fully correlated state 

^2c n (t 2 )\n) <g) \a n ) <g> \e n ). (30) 

n 

The final step is to argue that the environment states |e n ) are by definition inaccessible 
or uninteresting, so one should trace over them to obtain the reduced density matrix for 
system-plus-apparatus. Doing this we obtain 

=^2\c n {t 2 )\ 2 \n; a n )(n; a n |, (31) 

n 

where we have taken the environment states to be orthonormal, and we have defined 
\n; a n ) = \n) ® \a n ). This reduced density matrix describes a mixture of states which 
have classically correlated system-apparatus variables. The measurement is complete, with 
all quantal correlations between system and apparatus having diffused into the environment, 
thus effectively becoming lost. 

We will now connect this schematic to our case by examining a microscopic energy- 
changing collision between a system particle and a spin-0 bath boson 0. We previously called 
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the Hilbert space spanned by the horizontal and vertical states as the "system" . It is now 
useful to change perspective. Let us redefine "system" to mean the Hilbert space spanned 
by just the horizontal eigenstates (left or right), with the "apparatus" being spanned by the 
vertical states (single well energy eigenstates labeled by E). The environment is spanned by 
the energy states |e) of <fi. 
The initial state is 

[ Cl (0)\L,E) + c 2 (0)\R,E)]®\e), (32) 
the left and right eigenstates being 

\L,E) 



\E + u) + \E-u) 



V2 
) - |- 
V2 

respectively, where E±u label the nearly degenerate double-well states collectively indexed 



\E + u) - \E -u) 
\R,E) = '-, (33) 



by the vertical designation E. Free unitary evolution takes the inital state in Eq. fl32|) to 

e-^ E+ ^[ Cl (t)\L, E)+c 2 (t)\R, E)] ® \e) (34) 

at time t, where Ci^{t) = Ci t2 {0) cosut — ^1(0) sin ut. We will ignore the unimportant 
overall phase factor in the following. Suppose a collision with a bath boson occurs at the 
time t (we are assuming that the interaction time is very short compared to the mean free 
path). For simplicity, consider the extreme case b = 1 so that only L states interact with (/>, 
and work in one spatial dimension. The latter idealisation, which is inessential, leads to the 
final energies being completely determined by kinetic energy and momentum conservation 
from the initial energies E and e. Immediately after the collision, the state of the system is 

cx(t)\L, E') ® \e'} + c 2 (t)\R , E) ® |e), (35) 

where E' 7^ E and e' 7^ e.[] This is already in the fully entangled form of Eq. (|30|) ; there 



7 In addition, the oscillation frequency has changed from u> to cu' for the left eigenstate. 
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is no distinct pre-measurement step. We now trace over the energy states to obtain the 
reduced density matrix 

| Cl (t)| 2 |L, E')(L , E'\ + \c 2 (t)\ 2 \R, E)(R, E\. (36) 

In the sense made precise by the comparison of the above simple argument with the general 
schema, the measurement of the horizontal variable is complete. 
On the other hand, if b = 0, then Eq.(|35[) is replaced by 

[ci(t)\L, E') + c 2 (t)\R, E'}] ® \e'}, (37) 

whose reduced density matrix is that for the pure state C\(t)\L , E') + c 2 (t)\R, E'). There is 
no entanglement with the environment. The horizontal quantity remains unmeasured, even 
though there has been an energy changing collision with a bath boson. 

It is interesting that it makes sense to consider the vertical states as spanning an appara- 
tus which measures the horizontal states. We began by considering a three-part division of 
the world (horizontal, vertical, bath) for other reasons, but apparently this same partition 
provides a framework for interpreting aspects of the 67^0 dynamics as a canonical measure- 
ment process. Indeed, our master equation incorporating both horizontal and vertical states 
then precisely specifies the dynamics of what can be interpreted as a system-plus-apparatus 
hybrid. 

We can sharpen this point. Suppose we put some particles into one of the energy eigen- 
state of an isolated well. They oscillate between the left and right states at a given frequency. 
At time tbath (bath time) suddenly immerse it in a 6 = 1, Q — >oo bath. If the particles' 
phases happen to place them entirely in the left state at this time, then the system quickly 
becomes thermalised in energy and frozen in the left eigenstate. Now remove the system 
from the bath. If a thermalised distribution of oscillating particles is seen, then we know 
that the horizontal state was measured to be "left". On the other hand, if the phases at 
tbath happen to place all of the particles in the right eigenstate, then the system will remain 
unchanged in energy but frozen in the right eigenstate. Upon removal from the bath, we see 
a non-thermalised distribution and hence we know the measurement yielded "right" . 
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V. CONCLUSION 



We have analysed behaviors that can be exhibited by a quantum system that is coupled 
to a thermal bath. The system states were taken to have both a horizontal and a vertical 
structure. The horizontal states are discrete and the system is generally in a superposition 
of them. Familiar examples are the left and right sides of a symmetrical double well, the 
case we focussed on, and neutrino flavor. The vertical states are labeled by the energy 
eigenvalues. A "dynamical superselection rule" induced by decohering interactions with the 
bath quickly eliminates superpositions in the vertical direction. However, the same system- 
bath coupling induces system transitions between the vertical states. The limit where the 
system-bath interaction rates were much larger than the free oscillation rates of the system 
was our particular concern. 

Two qualitatively different, but complementary, limiting cases were analysed. The frozen 
or quantum Zeno case arises when the system-bath interactions distinguish between the 
horizontal states (side-sensing). The synchronized or motionally narrowed case arises, by 
contrast, when the system-bath interactions are blind to the horizontal structure (side-blind). 
The latter can only be observed when the vertical structure is present, whereas the former 
obtains even when the system is vertically trivial. Two extensions were made: to a double- 
well containing fermions, where the effect of Fermi statistics was fully incorporated into 
the master equation; and to potentials with a time-dependent bias. In the synchronization 
limit, the multi-fermion system presents a behavior which gives the illusion that all of the 
fermions are in the same state, evolving coherently. When there is a time-dependent bias, 
the synchronization limit can lead to an increase in the efficiency of transformation from 
one horizontal state into the other. 

We then examined the time evolution of bipartite entanglement-of-formation in the face 
of side-blind system-bath interactions. Entanglement was observed to at first decrease with 
increasing system-bath coupling strength. However, when the synchronization limit was 
reached, entanglement was found to be preserved. 
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The connection between freezing/QZE on the one hand, and synchronization/motional- 
narrowing on the other, and measurement was clarified. By adopting the emerging "system 
+ apparatus + environment" paradigm for quantal measurement, we were able to show that 
microscopic energy-changing interactions between a system particle and a bath boson can 
be naturally interpreted as measurement events in the horizontal observable. The tripartite 
subdivision "horizontal + vertical + bath" was mapped on to "system + apparatus + 
environment" in the process. 

In closing, we would like to note that the synchronization or motional narrowing limit 
should be useful in applications such as quantum computation which exploit quantal coher- 
ence and entanglement. The experimental or engineering challenge would be to construct 
devices whose active or horizontal system elements are deliberately placed in a noisy envi- 
ronment, but one that couples to all horizontal states equally strongly. 
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APPENDIX A: DERIVATION OF THE MASTER EQUATIONS 

We provide a derivation of the master equations appropriate for the context of this paper. 
We focus on the aspects that give rise to the placing of the horizontal matrices in Eq.(f|) 
and the detailed form of the right hand side of Eq. ([TDD . The methods are closely related to 



those given in |l7j but we employ no Markov approximation. We introduce an interaction 



picture based on the above division of the Hamiltonian. The system-bath coupling in this 
picture is denoted by Hj(t), with vertical matrix elements, 

The leading terms in the system+bath density matrix, p} +h , in this interaction picture are 
generated by the integral equation, 

p s + b (i) = p s+h (t = 0) 

rt rh 



I dt x r dhiip^ih),^^)],^^)}. (A2) 

JO JO L J 



In the perturbation expansion of Eq.( |A2|) we find a piece with an additional power of the 
quantities dj k [Et, EX -1 ] in each higher order. This provides our definition of leading terms 
and the explanation for how a tiny c| fc can lead to appreciable effects over long periods of 
time. These terms are generated by the terms in which each successive new pair of if/s 
in the iteration solution puts the state back to the same vertical (and bath) state. These 
observations provide the basis for having omitted, in writing Eq.( [A2[ ), all odd terms in Hj, 
and they also dictate the form of the density matrix to be used in solving the equation, 
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(A3) 



where, 



bath 



l-e 



exp 



{^2ala n uj n )/T 



(A4) 



and where, 



p I (E ] ,t)=e^ t p(E J ,t)e-^ 



i\(Ej)t 



(A5) 



We substitute Eq. (|A3|) into Eq. (|A2|) , take the bath trace, and use leading order expressions 
for each term in the double commutator, obtaining, for example, 



x J2 r ( E k, E^e^^Xe-^^ 1 

k 

x Pl (E k My KEk)t Ke~ tKE3)tl ■ 



8{t x - 1 2 ) 



(A6) 



Similarly, we have 



Tr 



bath 



{jlHj^HjihW b (h)\j) = 5(h - t 2 ) 



x e 



k 



(A7) 



The key to separation of the leading terms in the above is performing the sum over the modes 
of the scalar field first. After making the transition to the continuum, we make replacements 
of the form, 



-> 27rf[(E k -E j )]5(tx-h)9(E k -E j ) 



(A8) 



where neither the introduced infrared cut-off, w min nor the residual terms contribute to 
leading order, as defined above. 

Using Eqs.( |A6| ), ( |A7| ) and their counterparts for the two other orderings in Eq.(6), doing 
the t 2 integral (with the delta function symmetrically smeared, as can be justified by more 
accurate integrations) and differentiating with respect to t gives, 
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I Pl {Ej, t)=Y, e^\e-^ E ^ Pl {E k , t) 
x e iA(£fc)t Ce- iA( ^ )t r(E fc ,E i ) 

_l^ e iX(E j )t C 2 e -iX(E j )t pi ( E ^ t > ) 

+ Pl (E J ,t)e lX ^ t C 2 e- lX ^ t )J2nE 3 ,E k ). (A9) 
Using Eq. (|X5|) we now regain Eq.([|). 
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